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Abstract 

This paper presents the design of a multi-spacecraft system for the deflec- 
tion of asteroids. Each spacecraft is equipped with a fibre laser and a solar 
concentrator. The laser induces the sublimation of a portion of the surface of 
the asteroid, and the resultant jet of gas and debris thrusts the asteroid off its 
natural course. The main idea is to have a formation of spacecraft flying in the 
proximity of the asteroid with all the spacecraft beaming to the same location 
to achieve the required deflection thrust. The paper presents the design of the 
formation orbits and the multi-objective optimisation of the formation in order 
to minimise the total mass in space and maximise the deflection of the aster- 
oid. The paper demonstrates how significant deflections can be obtained with 
relatively small sized, easy-to-control spacecraft. 
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1. Introduction 

Near Earth Objects (NEO), the majority of which are asteroids, are defined 
as minor celestial objects with a perihelion less than 1.3 AU and an aphelion 
greater than 0.983 AU. A subclass of these, deemed Potentially Hazardous Aster- 
oids (PHA) , are defined as those with a Minimum Orbital Intersection Distance 
(MOID) from the Earth's orbit less than or equal to 0.05 AU and a diameter 
larger than 150 m (equivalent to an absolute magnitude of 22.0 or less). As of 
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March 2012, 8758 NEO's have been detected (IAU Minor Planet Centre, 2012); 
of those, 840 are estimated to have an effective diameter larger than 1 km*, 
and 1298 are categorised as potentially hazardous. Impacts from asteroids over 
1 km in diameter arc expected to release over 10 5 megatons of energy with 
global consequences for our planet (Stokes et al., 2003), while those with an 
average diameter of 100 m can are expected to release over 10 2 megatons of en- 
ergy potentially causing significant tsunamis and/or land destruction of a large 
city (Toon et al., 1997). It is estimated that there are between 30,000-300,000 
NEO's with diameters around 100 m, meaning a large number of NEO's are still 
undetected. 

A quantitative comparison of the various options for NEO deflection was 
conducted by Colombo et al. (2006); Sanchez Cuartielles et al. (2007). Examin- 
ing the results of the comparison, one of the more interesting methods employed 
solar sublimation to actively deviate the orbit of the asteroid. The original con- 
cept, initially evaluated by Mclosh et al. (1994), and later assessed by Kahlc 
et al. (2006), envisioned a single large reflector; this idea was expanded to a 
formation of spacecraft orbiting in the vicinity of the NEO, each equipped with 
a smaller concentrator assembly capable of focusing the solar power at a dis- 
tance around 1 km and greater (Maddock et al., 2007). This concept addressed 
the proper placement of the concentrators in close proximity to the asteroid 
while avoiding the plume impingement and provided a redundant and scalable 
solution. However, the contamination of the optics still posed a significant lim- 
itation as demonstrated by Vasile and Maddock (2010). In the same paper, the 
authors demonstrated that the combined effect of solar pressure and enhanced 
Yarkovsky effect could lead to miss (or deflection) distances of a few hundred 
to a thousand kilometres over eight years of deflection time. However, this de- 
flection is orders of magnitude lower that the one achievable with a prolonged 
sublimation of the surface. 

A possible solution is to use a collimating device that would allow for larger 
operational distances and protection of the optics. This paper presents an as- 
teroid deflection method based on a formation of spacecraft each equipped with 
solar pumped lasers. The use of lasers has already proposed by several authors, 
although always in conjunction with a nuclear power source (Phipps, 1992, 1997; 
Park and Mazanek, 2005). Extensive studies on the dynamics of the deflection 
with high power lasers were proposed by Park and Mazanek (2005) envisaging a 
single spacecraft with a MW laser. This paper proposes a different solution with 
a formation of smaller spacecraft, each supporting a kW laser system indirectly 
pumped by the Sun. 

The paper starts with a simple sublimation model that is used to compute 
the deflection force. The orbits of the spacecraft formation are then designed 
by solving a multi-objective optimisation problem that yields an optimal com- 
promise between distance from the target and impingement with the plume of 



tAn asteroid with an effective diameter equal to or greater than 1 km is defined here to be 
any NEA with an absolute brightness or magnitude H < 17.75, as per Stuart (2003). 
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debris. A Lyapunov controller is proposed to maintain the spacecraft in for- 
mation along the desired proximal orbit. A second multi-objective optimisation 
problem is then solved to compute a different type of controlled formation orbits 
in which the shape of the orbit is predefined. Finally, the number and size of 
the spacecraft is optimised to yield the maximum possible deflection. 

2. Deflection Model 

The orbital properties of Near Earth Asteroids (NEA) can be grouped into 
four general categories based on the semi-major axis a of the orbit, radius of 
apoapsis r a and/or radius of periapsis r p , described as follows (NASA Near 
Earth Object program, 2012): 

Atens Earth-crossing asteroids with semi-major axes smaller than Earth (named 
after asteroid 2062 Aten), where a < 1 AU, r a > 0.983 AU. 

Apollos Earth-crossing asteroids with semi-major axes larger than Earth (named 
after asteroid 1862 Apollo), where a > 1 AU, r p < 1.0167 AU. 

Amors Earth-approaching asteroids with orbits exterior to Earth's but interior 
to Mars (named after asteroid 1221 Amor), where a > 1 AU, 1.0167 AU 
< r p < 1.3 AU. 

Atiras Near Earth Asteroids whose orbits are contained entirely with the orbit 
of the Earth (named after asteroid 163693 Atira), where a < 1 AU, r p < 
0.983 AU. 

Apollo is the largest class (approximately 4100 NEA's) followed by Amors (ap- 
proximately 3400 NEA's). The asteroid Apophis 99942, part of the Apollos 
class, is taken as a test case with a relatively low aphelion such that enough 
solar power can be harvested. 

While circular, or near-circular, orbits offer a more constant level of solar ra- 
diation, as suggested by Vasile (2008); Vasile and Maddock (2010) if the mirrors 
have variable optics, i.e., the focal point can be changed, a constant power den- 
sity can be achieved for asteroids on elliptical orbits or when the level of solar 
power available is low. In terms of altering an orbit, thrusting at the perihelion 
of elliptical orbits maximises the change in semi-major axis (and therefore the 
miss distance). Even if the level of solar radiation available is not sufficient 
to induce sublimation at aphelion, a deflection can still be achieved, as will be 
demonstrated in this paper. 

The other benefit of basing the test case on the Apophis asteroid is its 
popularity in scientific literature due to the initial, relatively high impact level 
(2.7% chance of impacting the Earth in 2029) it was given when it was first 
observed in 2004. While further tracking data has reduced the threat level, 
ruling out the possibility of an impact in 2029 but leaving a non-zero impact 
probability for the 2036 and 2037 encounters, the asteroid Apophis remains a 
popular reference example. Note, in fact, that although Apophis is not necessary 
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a typical case, the interest here is to examine the effectiveness of a fractionated 
laser ablation system applied to the deflection of an S-type asteroid, belonging 
to a given size range, on a moderately eccentric orbit (although also an extension 
to highly eccentric orbits will be demonstrated) of which Apophis is an example. 

Table 1 gives the orbital and physical data of the asteroid used in this study. 
The asteroid shape was assumed to be tri-axial ellipsoidal, 

at = V2d A b t = d A c e = — ^ (1) 

v/2 

where at > be > eg are the three radii along the three orthogonal axes and 
d A is the estimated average diameter based on the observed magnitude, given 
in Table 1. S-type asteroids, as used here for the test case, are moderately 
bright with an albedo from 0.10 to 0.22. By comparison, C-type asteroids are 
extremely dark with albedos typically in the range of 0.03 to 0.10. According to 
Delbo et al. (2007), the geometric albedo of Apophis is 0.33 however the value 
used here of 0.2 was chosen to give a more general test case. 



Table 1: Orbital and physical properties of test asteroid. 



Element 




Measured Value 


Semi- major axis 


a A 


0.9224 AU 


Eccentricity 


e Ao 


0.1912 


Inclination 


'a„ 


3.3312 deg 


RAAN 




204.4428 deg 


Argument of pcriapsis 


W A 


126.4002 deg 


Period 


T Aq 


323.5969 days 


Mean motion 


n Ao 


1.2876 xl0~ 5 deg/s 


Mass 


m Ao 


2.7xl0 10 kg 


Gravitational constant 


fJ-A 


1.801599xl0~ 9 km 3 /s 2 


Physical dimensions 


a,£,be,C£ 


191 m, 135 m, 95 m 


Rotational velocity 


W A 


3.3xl0~ 3 deg/s 


Albedo 




0.2 



The minimum orbital intersection distance (MOID) is the separation dis- 
tance at the closest point between two orbits, e.g., Apophis and the Earth. 
The deviation distance is defined here as the difference in position between the 
original, undeviated orbit k Ao and the deviated orbit k Adeij at £moid (Colombo 
et al., 2009b) (see Fig. 1). Figure 2 illustrates the reference frames used here, 
where 0{i, j, k} is the inertial heliocentric reference frame, and A{x, y, z} is the 
relative, rotating Hill reference frame (radial x, transverse y and out-of-planc z 
directions), centred on the asteroid. 

Non-linear equations were used for determining the asteroid deviation vec- 
tor Ardev — r A de „ — Tao as a function of the ephemeris in the Hill reference 
frame A, as derived by Maddock and Vasile (2008), where Ak = k Adcu — k Af) = 
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Thrust Arc 



Figure 1: Definition of deviation distance at Figure 2: Definition of the reference frames, 
the MOID. including the rotating Hill frame A centred 

on the asteroid. 



[Aa, Ae, Ai, Af2, Auj, AM] T giving the difference in Kcplerian parameters be- 
tween the undeviated and deviated orbits. 

The change in the orbital parameters is calculated by numerically integrating 
the Gauss planetary equations (see e.g., Battin, 1999) using a thrust vector 
Udev = [ih u n Uh] T in the tangential, normal and out-of-plane (or direction of 
angular momentum h) reference frame, induced by the deflection method: 



Within this study, the deflection action is assumed to be aligned with the he- 
liocentric velocity of the asteroid, therefore u n = and Uh = 0. Other authors 
have studied the optimal direction of the deflection action in the case of laser 
ablation (Yoo et al., 2009), however, the main interest of this paper is in the 
system sizing in relation to the achievable deviation. 

Colombo et al. (2009b) determined that the change in angular location, in 
this case given by the mean anomaly M, calculated at the MOID is, 

AM = dt + n Ag (t - t MOlD ) + n Ai (t MOID - U) (3) 

where Ua is the mean motion of the undeflected asteroid, is the mean 
motion of the asteroid at the end of the deflection action, to is the beginning of 
the deflection action and ti is the end of the deflection action. 

The non-linear proximal motion equations in Vasile and Maddock (2010) 
together with Eq. (3) and the Gauss planetary equations give the variation of 
the orbit of the asteroid at the time of the MOID. Vasile and Colombo (2008) 
showed that an estimation of the minimum orbit interception distance can be 
computed by projecting the variation of the orbit at the expected impact time 
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onto the 6-plane of the Earth at the time of the MOID, i.e., computing the 
variation of the impact parameter b. Hence, in the test section the variation of 
the impact parameter will be used as a measure of the achievable deflection. 

The thrust produced by the deflection method is computed assuming that 
the lasers are not pulsed but continuous wave and that the energy density is 
sufficient only to turn the matter into gas (vapour regime) but not to produce 
plasma (Phipps, 2010). The level of momentum coupling that can be achieved 
with this model is lower than what can be found in other studies (see e.g., 
Phipps, 2010). A further assumption is that the asteroid is absorbing part of 
the incoming energy without changing its temperature thus providing a constant 
sink for heat transmission; this might not be the case for small asteroids. 

Under these assumptions, the rate of the expelled surface matter is defined 
as (Sanchez Cuartielles et al., 2009), 

dfTlexp ~ f f 



dt 



rVmax flout Y 

2n sc v rot / / — {Pi„ - Qrad - Qcond) dt dy (4) 
Jyo Jt in H 

where [ti n ,t ou t] is the duration for which a point is illuminated, [yo,Vmax] are 
the vertical limits of the illuminated surface area (i.e. orthogonal to the direction 
of rotation of the asteroid), H is the enthalpy of sublimation, v rot is the linear 
velocity of a point as it travels horizontally (i.e., orthogonal to y) through the 
illuminated spot area and n sc is the number of spacecraft in the formation. 
The input power per unit area due to the solar concentrators is given by, 

Pin = VsysC r (l - <; A )S C~f^j ( 5 ) 

where <^ A = 0.2 is the albedo, Sq = 1367 W/m 2 is the solar flux at 1 AU, scaled 
to the Sun-asteroid distance r A , rj sys is the system efficiency, and C r is the 
concentration ratio (the ratio between the power density from the Sun on the 
mirror surface, and that of the illuminated spot area on the asteroid). 

The heat loss due to black-body radiation and the conduction loss are de- 
fined, respectively, as, 

Qrad = crebbT 4 (6) 



Qcond — {T su bl — ^o)y A (7) 

where a is the Stefan-Boltzmann constant, ej,b is the black body emissivity, 
T is the temperature and c A , Pa and k A are, respectively, the heat capacity, 
density and thermal conductivity of the asteroid. For the asteroid Apophis, c A = 
750 J/kg-K based on the average value for silicate materials, k A = 2 W/K/m 
and p A = 2600 kg/m 3 (Remo, 1994). The sublimation temperature assumed is 
that for forsterites (Wang et al., 1999), T suU = 1800 K, with T set to 278 K. 

The induced acceleration due to the sublimation process can then be deter- 
mined by (Sanchez Cuartielles et al., 2009), 

Avm exp „ 

u su6 = ^v A (8) 
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where ttia is the mass of the asteroid at a generic instant of time, v A is the 
direction of the velocity vector of the NEO, A ~ (-) is the scattering factor, v 
is the average velocity of the debris particles according to Maxwell's distribution 
of an ideal gas: 



where k B is the Boltzmann constant, and MM S ,si0 4 is the molecular mass of 
fosterite. 

The scattering factor A is computed as the average of all possible thrust di- 
rections assuming that the thrust can point randomly at any angle a T between 
and 7r, therefore A = i J Q cosa T do; T (Sanchez Cuarticlles ct al., 2009). Some 
preliminary experiments (Gibbings et al., 2011) demonstrate that the plume 
is progressively focusing inwards for rocky type of asteroids, while for highly 
porous asteroids the plume tends to remain unfocused; hence assuming an uni- 
form distribution of the thrust pointing direction over an angle of 180° is a 
conservative choice. The remaining mass of the asteroid m A is calculated by 
numerically integrating Eq. (4). 

2.1. Contamination Model 

The contamination of the mirror surfaces due to the debris plume is modeled 
based on the work by Kahle et al. (2006). Their study made a number of initial 
assumptions regarding the expansion of the plume and sublimation process. 
The first assumption holds that the sublimation process is comparable to the 
generation of tails in comets. The asteroid is assumed to contain a reservoir 
of material underneath the surface, with the gas expanding outwards through 
a throat into vacuum. Preliminary experimental results have shown that this 
assumption, as with others in this section, are potentially overly pessimistic and 
may not be valid for every type of asteroid. However, altering these assumptions 
does not change the fundamental results in this paper, therefore it was decided 
to remain consistent with the existing literature and defer any further analysis 
on the validity of these assumptions for future work. 

The second assumption is that the plume expansion is similar to the expan- 
sion of gas of a rocket engine outside the nozzle. The density of the expelled 
gas p exp is computed analytically, 

Pe Xp (5r s/sc , <p) = 3c L dsP °! , ) 2 (cos e) 2 ^" 1 ) (10) 

v A spot \ 2dr s/sc + d spot J 

where d spo t is the diameter of the spot area, 5r s / sc is the distance from the spot 
on the surface of the asteroid and the spacecraft, and 6 = irtp/2(p max where ip is 
the angle between the spot-spacecraft vector and the y-axis of the Hill reference 
frame. The jet constant j c was set to 0.345, the maximum expansion angle 
Vmax = 130.45°, and adiabatic index n = 1.4 based on the values for diatomic 
particles (Legge and Boettcher, 1982). 
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Note that this density model is in contradiction with the assumption of a 
uniform scattering over a hemisphere and, in fact, suggests a much more focused 
plume. From ongoing experiments (Gibbings et al., 2011), the plume appears 
to more closely match the density distribution given in Eq. (10) rather than 
a uniform distribution; nevertheless, in the analysis in this paper the most 
conservative choice was selected for the scattering factor in order to account for 
possible unmodclcd performance degradation components. 

The position vector <5r s / sc from the spot to the spacecraft is defined as: 



x — ri sinw A t cos(— w A t — 9 V J + ri cosw A tsm(—w A t — 9 Va ) 
y — ri cosw A tcos(—w A t — 9 V J — rg sin w A t sin(— w A t — 9 Va ) 

z 



(11) 



where the radius of the ellipse is given by, 



^ = I 2 =f ( } 

\J (be cos(-w A t - 9 Va )) + (a f sin(~w A t - 9 V J) 

and, with reference to Fig. 2, the position of the spacecraft with respect to the 
centre of the asteroid is 6r = [x, y, z] T . We assume here that the asteroid 
is spinning around the z axis with a rotational velocity w A . The direction of 
the velocity of the asteroid in the heliocentric reference frame projected onto 
the Hill reference frame A is 6 Va . In other words, in order to have a deflection 
thrust aligned with the velocity of the asteroid, the spot is assumed to be at an 
elevation angle over the y-axis equal to 8 Vl . 

The third assumption made is that all the particles impacting the surface of 
the mirror condense and stick to the surface. The exhaust velocity is constant, 
therefore the thrust depends only on the mass flow. A higher thrust results 
in a higher mass flow and thus in a faster contamination. This is a rather 
conservative assumption. The actual contamination level depends on the type 
of deposited material and the temperature of the optical surfaces. Following the 
approach used to compute the contamination of surfaces due to out-gassing, a 
view factor ip v f was added equal to the angle between the normal to the mirror 
and the incident flow of gas. The resulting variation of the thickness of the 
material condensing on the mirror can be computed by, 

^nd = ^Pexp cQsijvf (13) 

"™ Player 

The average debris velocity v is multiplied by a factor of 2 to account for the ex- 
pansion of the gas in a vacuum. The layer density p la er was set to 1 g/cm 3 . The 
power density on the asteroid surface is decreased based on the contamination 
of the mirrors. 

A degradation factor r is applied to the power beamed to the asteroid surface, 
based on the Lambert-Beer-Bouguer law (Kahle et al., 2006), 

t = e - 2vh ^ d (14) 



where v = 10 4 /cm is the absorption coefficient for forsterite. Note that the 
values of v and p laycr are based on the assumption that the deposited material 
is dense and absorbs the light over the whole spectrum. This is again a rather 
conservative assumption; experiments have shown that while it appears to be 
valid for some silicates such as forsterite, this assumption may not hold true for 
all materials. As mentioned previously, further experimentation and analysis 
are underway, and will be the topic of future publications. 

Eq. (13) is numerically integrated, along with the Gauss equations, for the 
period of the mission. 

2.1.1. Tugging Effect 

The spacecraft will fly in formation with the asteroid at a distance 5r, thus 
exerting a gravitational pull on it (Gong et al., 2009). The tugging acceleration 
utug is given by: 

utug = -n sc -^-6r (15) 

where G is the universal gravity constant and m sc is the mass of a spacecraft. 
The sum of u tug and u su b forms the total deflection acceleration \idev The 
acceleration \idev is used with Gauss planetary equations in order to determine 
the change in the NEO orbit. 

2.2. The Laser System 

Lasers work on the general premise of exciting electrons by stimulating them 
with the addition of photons (or quantum energy), which temporarily boost 
them up to a higher energy state. This stimulation continues until a population 
inversion exists, where there are more electrons at a higher energy state, e.g., 
Ei than at the lower (or original) state, e.g. Eq. The release of photons when 
the electrons drop back to their original base state produce an emission that, 
generally, has the same spectral properties of the stimulating radiation, and 
is therefore highly coherent. The energy that is not released as part of the 
output emission, is instead released as heat. This means that the laser must be 
continually cooled, which in space means large radiators. 

In this paper two general methods of powering the laser are considered and 
defined as: direct pumping, where the energy is directly used to excite the laser, 
and indirect pumping, where an intermediate step is used to first convert the 
energy, e.g., solar radiation, into electricity. 

Indirect solar-pumped lasers convert the solar energy first into electricity, 
which is then used to power the laser. Photovoltaic cells are an obvious choice 
for space applications. The drawback, of course, is the addition of an electri- 
cal power generator meaning added mass, size and power requirements. Direct 
solar-pumped lasers, by comparison, do precisely what the name suggests: the 
laser is directly energised using solar radiation. Due to the mismatch between 
the wide-band emissions of the Sun with the narrow absorption bands of lasers, 
the loss of available solar power is currently rather high. For example, the over- 
lap between a Nd:YAG (neodymium-doped yttrium aluminium garnet) crystal 



absorption spectrum and the solar radiation spectrum is around 0.14 (Weksler 
and Shwartz, 1988). 

One option is to use high efficiency solar arrays in conjunction with a solid 
state laser. Solid state lasers pumped with electric power can currently reach 
60% efficiency. If the solar arrays have an efficiency of 30%, then the system 
would have an overall efficiency of 18%. If a pumped laser is used, then the 
focal point can be close to the primary mirror and a high concentration factor 
can be obtained with a relatively small mirror. For example, if the mirror has 
an area of 314 m 2 (equivalent to a 20 m diameter circular mirror), then the 
collected power at 1 AU is 429.5 kW. The solar array plus laser system converts 
only 18% of this power, therefore only 77.3 kW is beamed to the surface to the 
asteroid, while the rest needs to be dissipated. 

In a paper presented in 1994, Landis discussed the use of a directly solar 
pumped laser based on semiconductor technology. According to Landis, the 
expected efficiency of directly pumped semiconductor laser would depend on 
the same efficiency losses of a solar cell, therefore Landis was expecting a lasing 
efficiency (output/input power ratio) of 35%. Such an efficiency would be one 
order of magnitude higher than the best Nd:YAG laser system, which is expected 
to reach 6% of overall efficiency. 

Direct solar pumping would represent an interesting solution in terms of 
complexity of the overall system. In fact no cooling system for the photovoltaic 
conversion and no power transmission would be required. On the other hand 
the Technology Readiness Level (TRL) of both solar cells and semiconductor 
lasers is far higher than the one of a directly pumped laser and an indirectly 
pumped laser can be expected to be operational much sooner. 

Recent electrically pumped semiconductor laser have proven over 73% wall- 
plug efficiency (Crump ct al., 2005; Stickley et al., 2005; nLIGHT. 2006; Peters 
et al., 2007) with a target efficiency of 80%. Research on fibres coupled with 
clusters of diodes have demonstrated slope efficiencies of up to 83% (Jeong et al., 
2003, 2004). A substantial increase in cells efficiency has also to be expected. In 
particular, in order to achieve a 35% efficiency in direct pumping, semiconductor 
technology should allow the absorbtion of the solar spectrum over a wide range of 
frequencies. A high efficiency of a directly pumped laser is therefore expected to 
correspond to a high efficiency of solar cells. An increase of solar cell efficiency 
up to 50% (Luque et al., 2004) is reasonable, allowing an indirect pumping 
system to have a comparable efficiency to a 35% direct pumping system. 

In the following the assumption is that the overall system efficiency r) sys is 
about 22.7%, with a 45% efficiency of the cells, a 90% efficient reflectors, a 85% 
efficiency of the power transmission and regulation line and a 66% efficiency of 
the laser (given by the product of the target 80% for the laser diode and the 
achieved 83% slope efficiency of the fibres). A second option with a 60% laser 
efficiency and 40% cell efficiency is also considered. 
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3. Formation Design 

One idea for the orbital design is to have the spacecraft flying in forma- 
tion with the asteroid, orbiting in tandem around the Sun (see Fig. 11). The 
spacecraft have to maintain their relative position with respect to the asteroid 
in order to keep the required power density on the same spot on the surface 
of the asteroid (note that the surface of the asteroid is moving under the spot 
light of the laser). Therefore, the formation orbits have to be periodic and in 
close proximity with relatively low excursion in the relative distance from the 
asteroid. On the other hand the spacecraft should minimise any impingement 
with the plume of debris and gas coming from the sublimation of the surface 
material. 

In order to design the desired formation orbits, one can start by considering 
the local Hill reference frame A{x,y,z} in Fig. 2 and the associated linearised 
version of proximal motion equations (Schaub and Junkins, 2003) used in the 
calculation of the asteroid deviation vector: 

x[y) = 6M — aACOsvbe (16a) 

V 

t* f sin v 
y{y) = — (1 + eA cos v) 2 5M + r K 5uj H ^ — (2 + e A cosv)5e + r A cos iaSQ 



where rj = yl — e\, 6 A = v + oja, v is the true anomaly, a a, e-A, iA, 
are respectively the semi-major axis, eccentricity, inclination and argument of 
the perihelion of the orbit of the asteroid at a generic moment in time and 
5k. = [5a, 5e, Si, <5f2, Soj, 5M] T are the variations of the orbital elements, with 
the imposed conditions 5r <C r A , and 5a = in order to have periodic motion. 
These equations are a first approximation of the motion of the spacecraft and 
do not take into account the gravity field of the asteroid or solar pressure. 

If the optimal thrust direction that maximises the deviation is along the 
unperturbed velocity vector of the asteroid (Colombo et al., 2009b), then the 
exhaust gases will flow along the direction of the velocity of the asteroid pro- 
jected in the Hill reference frame. Therefore, the position vector in the radial, 
transversal and out-of-planc reference frame was projected onto the tangential, 
normal, out-of-plane reference frame to give 5r tn h = [xtnh,ytnh, z tn h] T ■ Then, 
the size of the formation orbits projected in the Xtnh-ztnh plane was maximised. 
All the requirements on the formation orbits can be formulated in mathematical 
terms as a multi-objective optimisation problem with two objective functions, 



z{y) 



r A (sin 6a5i — cos Q a sin Ia5£1) 



(16b) 
(16c) 



min max J\ = 5r 

5kG-D v 



(17a) 




(17b) 
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Tabic 2: Boundaries on the formation orbital parameters 





Se 


Si 


SQ 


Suj 


5M 




(io- 7 ) 


(10~ 7 rad) 


(10~ 7 rad) 


(10~ 7 rad) 


(10- J rad) 


Lower bound 


-0.01 


-0.1 


-0.9 


-1.5 


-0.1 


Upper bound 





0.1 


0.9 


1.5 


0.5 



subject to the inequality constraint, 

C ineq = min \y{v)\ - y U u > (18) 

where y hlM is a minimum distance along the y-axis, and D is the search space 
for the solution vector (5k. Table 2 defines the boundaries imposed on D. The 
boundary values were obtained by progressively increasing each of the bound- 
aries from to the value in the table, looking at the value of the maximum 
distance from the asteroid. Larger boundaries produce solutions with a better 
(lower) performance index J2 but a higher performance index J\. 

Equations (17)-(18) were optimised using a memetic multi-objective opti- 
miser MACS (Multiagent Collaborative Search) (Vasile, 2005; Maddock and 
Vasile, 2008; Vasile, 2008). The optimisation led to the identification of two 
families of formation orbits belonging to two subsets of the search space D. Fig- 
ures 3, 4 and 5 show the two families in the parameter space for y UM = 1000 m. 
The solutions are almost perfectly symmetrically distributed about the 0-value 
of 5i, Sfl, while there is a bias towards the negative axis for Suj. Each family 
has been identified with the label — z or +z depending on whether the sign of 
the z coordinate is negative or positive at y = j/ LIM . Figure 6, instead, shows the 
Pareto fronts for y LIM = 500 m and y UM — 1000 m respectively. Note that in Fig. 
6, the Pareto fronts for the branches in Figs. 3, 4 and 5 appear superimposed 
and cannot be distinguished. Therefore, the two families can be considered 
equally locally Pareto optimal. 

Figure 7 shows the formation orbits in the A Hill frame. It can be noted that 
the two families are symmetric with respect to the x-y plane. In the remainder 
of the paper these orbits will be called natural formation orbits. 

3.1. Formation Dynamics and Control 

In order to maintain the orbits designed in the previous section, the space- 
craft need to be controlled. In the proximity of the asteroid, in a Hill rotating 
reference frame, the spacecraft are subject to the force due to solar pressure, 
the gravity of the asteroid, the gravity of the Sun, the centrifugal and Coriolis 
forces plus the forces induced by the impingement with the plume. An active 
control is therefore required to maintain the spacecraft flying in formation with 
the asteroid. 

Following the Jacobi ellipsoid model, the minor axis ci of the asteroid is 
aligned with vector of angular momentum, which corresponds to the z-axis of 
the asteroid Hill frame A (see Fig. 2). The gravity field of the asteroid is 
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Figure 3: a) Objective functions J\ and J2 versus Si, b) objective function J2 versus 8i 
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Figure 7: Formation orbits with minimum distance of 500 m. 



expressed as the sum of a spherical field plus a second-degree and second-order 
field (Hu and Scheeres, 2002; Rossi et al, 1999), 

C/20+22 - £^ (c 2Q (1 - | cos 2 7 ) + 3C 22 cos 2 7 cos 2a) (19) 

where 7 is the elevation over the x — y plane and the harmonic coefficients C20 
and C22 are a function of the semi-axes, 

C 20 = -^(2c 2 -a 2 -& 2 ) (20a) 
Oa = ^(a?-6?) (20b) 

and A is defined as, 

A = arctan ^— J + w A t 

Note that a different rotational state or shape would alter the time- varying grav- 
ity field that the spacecraft would experience. In a real scenario, the rotational 
state coupled with the shape of the asteroid would require an adaptive focusing 
of the laser beam as the distance of the spot from the source will change with 
time. Also the divergence of the plume will change as the laser carves a groove 
into the asteroid. However, within the assumptions in this paper a different 
rotational state and/or shape would not alter the main results. 

If one considers a Hill reference frame A centred in the barycentre of the 
asteroid (see Fig. 2), the motion of the spacecraft in the proximity of the asteroid 
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is given by: 



•2/ , ^ , «„. MsuN^A + a;) Ma _ , F Sj: (x,y,z) i dU 2 Q+22 



a;(;/) = + 2^y + v (r A + x) + vy = —7 x 



Sr 3 m sc dx 

(21a) 



. 2 Msun Ma 



»(!/) = -2VX U(V A + x)+V y-—y-—y + ^ + ^ 

(21b) 

Kv ) = -*Lz + El^llA + g%±gg (21c) 



<5r 3 m sc <9z 



with, 



(22) 

'A 

r A = ^ r A — + u dev ^ (23) 

r A 

The force term F s = [F Sx F s F Sr ] T is made of two contributions: light pressure 
from the emitted light from the laser F srp and the force due to the flow of gas 
and debris coming from the asteroid F p i ume . 

The force due to solar radiation F srp is defined as, 

F srp = 2r lsys A Ml ^ (—) cos 2 (i n steer + (1 - r)l)A Ul ^ (—) x (24) 

where c is the speed of light and A Ml is the cross section area of the primary 
mirror (see 11). The angle j3 is the half angle between the normal to the steering 
mirror n steer and the Sun-mirror vector (which is approximated by setting it 
equal to the Sun-asteroid vector). The second term in Eq. (24) takes into 
account a non-perfect reflection of the primary and secondary mirror. The 
reflectivity of the two mirrors is here assumed to be rj M = 0.90. The assumption 
is that the energy dissipated by the radiators is emitted uniformly in every 
direction and does not contribute to any change in the linear momentum of the 
spacecraft. 

If the flow rate per unit area at distance Sr spot is (2p exp (Sr spot , f)v) and all 
the particles stick to the surface of the mirror then the force F p i ume is: 

" plume 

= Ap exp (Sr spotl ip)v 2 A eq cos ip vf Sr s/sc (25) 

The flow rate depends on the power density and therefore on the distance from 
the Sun. The part of the spacecraft exposed to the plume and to the reflected 
light changes along the orbit and is irregular. In order to simplify the calcula- 
tions, the assumption adopted in this paper is that the total effect is equivalent 
to a flat surface with area A eq = A Ml and normal unit vector h eq such that the 
cross product (h eq , <5r s / sc ) = cos ip v f. 
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Given these equations, the resultant of all the forces acting on the spacecraft 
is not zero and in particular the difference between gravity and F s is a function 
of time. Therefore, an active control is required to maintain the position of the 
spacecraft with respect to the asteroid. 

If one assumes that solar pressure, the gravity of the asteroid, and the force 
due to the plume impingement are the main source of perturbation of the prox- 
imity motion of the spacecraft and that any non-spherical terms in the gravity 
field expansion result in only a small (second order) additional perturbation, 
then one can build a simple control law based on the Lyapunov control func- 
tion: 

V = -5v 2 + -K \{x - x ref f + (y- y re ff + {z- z ref ) 2 J (26) 

where Sr re f = [x re f, y re f, z re f] T are the coordinates of a point along the nominal 
formation orbit (in the Hill frame A). The assumption here is that the motion 
along the reference formation orbit is much slower than the control action, which 
is valid as the period of the spacecraft orbit is equal to the period of the asteroid 
(just under 1 year). Therefore, the spacecraft targets a set of static points along 
the formation orbit. Now if there exist a control u such that dV/dt < then one 
can maintain the mirror in the proximity of the reference point as the reference 
point moves along the reference formation orbit. A possible control is given by: 

.ELs r + E^ + l^l)-K(Sr-Sr nf )-c d 6 V (27) 
6r 3 m sc m sc J 

The total derivative of the function V is: 

<Sv T <Sv + K(Sr - Sr ref ) T 6v (28a) 



dV 
~dt 



Sr 3 m sc m sc \ Sr 3 
— K (Sr — 5r re f) — c^v^ + K(5r — 5r re f) T 6v 

= -c d 5v T 5v < (28c) 

where 6v — [x, y, z] T is the relative velocity of the spacecraft in the asteroid Hill 
reference frame A. 

The control in Eq. (27) can now be introduced into the full dynamic model 
in Eq. (21) to test the validity of the assumption that the light coming from the 
asteroid and aspherical gravity field are indeed small. The elastic coefficient K 
for both cases was chosen to be 10 -6 while the dissipative coefficient Cd was set 
to 1CT 5 . 

Figure 8 shows the maximum thrust level as a function of the maximum 
distance from the asteroid for a 20 m diameter mirror. Figure 9 shows the 
propcllant consumption as a function of the maximum distance from the asteroid 
for a 20 m diameter mirror. 
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Figure 8: Maximum thrust versus maximum distance for the formation orbits with minimum 
distance of 1000 m. 



3.2. Shaped Formation 

Although the natural formation orbits are designed to minimise the impinge- 
ment with the plume of gas and debris, none of them can avoid the plume com- 
pletely. In order to maximise the amount of solar power collected, the mirrors 
should be constantly pointing directly towards the Sun, hence in a direction 
perpendicular to the y-axis. By following one of the natural formation orbits, 
the spacecraft will rise above the z-y plane (i.e, in the +x direction) once per 
revolution around the Sun, thus directly exposing the reflector to the plume. 
According to the contamination model, every surface directly exposed to the 
plume builds up a layer of a contaminants. This is quite a strong assumption as 
all the impinging material is assumed to condense and only the surfaces in view 
of the plume are contaminated. We hold on to these assumptions in this paper, 
although some experimental work is underway to build a more realistic model 
(Gibbings et al., 2011). If one sticks to the assumptions of the contamination 
model, then one solution to mitigate the contamination would be to fly always 
below the plume of gas (i.e., — x direction, below the z-y plane). In order to 
make the spacecraft follow the desired proximal motion the following shape is 
assigned to the formation orbit: 

x[y) = x\ cos(zy) + x-2 sin(f) + X3 (29a) 
y(v) = yi cos(v) + y 2 sin(zy) + y 3 (29b) 
z(v) = z\ cos(v) + Z2 sm(v) (29c) 

By differentiating with respect to time and inserting Eq. (29) and their first and 
second derivatives into the dynamic equations in Eq. (21), one can compute the 
control profile and the corresponding propcllant consumption. The interest now 
is to design formation orbits that minimise the propcllant consumption required 
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Figure 9: Propellant consumption versus maximum distance for the formation orbits with 
minimum distance of 1000 m. 



by the control system to remain below the z-y plane and operate as close as 
possible to the asteroid to minimise pointing requirements. The problem can 
be formulated as follows: 



min Ji = MF C 

sex 

min max Jo = 5?' 
sex v 

minmax J3 = ||u| 
sex v 



subject to the inequality constraints: 



C\ = maxx(p) < 

V 

C-2 = maxy(V) < 



(30a) 
(30b) 
(30c) 

(31a) 
(31b) 



where the solution vector is s = [x\, X2, X3, y-y, j/21 2/3 , z\, z 2 ] T , and MF C is the 
propellant mass fraction for the control over one year of operations. The search 
space X is defined by the lower and upper bounds on the components of s, re- 
spectively s ; = [-1,-1, -1,-1, -1,-2,-1, -If and s„ = [1, 1, 0, 1, 1, 0, 1, If. 
Again MACS was used to solve the constrained problem in Eq. (30) and Eq. (31). 

The result of the multi-objective optimisation can be found in Fig. 10, and 
shows the propellant mass fraction versus the maximum thrust level versus the 
maximum distance to the asteroid for the case of 10 spacecraft, each carrying a 
20 m diameter mirror, over the first year of operations. 

As expected the level of thrust and control propellant mass fraction are 
monotonically increasing with the distance to the asteroid. However, even for 
close distances the annual propellant consumption and the thrust level are quite 
small, only a few milli-Newton of thrust are enough to maintain the orbit. 
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Figure 10: Pareto front of the shaped formation problem. 



4. Spacecraft and System Sizing 

The proposed configuration of each spacecraft is as follows: each spacecraft 
is made of a primary mirror that focuses the sunlight onto a secondary mirror 
that reflects the light onto a solar array seated behind of the primary mirror (see 
Fig. 11). The electric power coming from the solar array pumps a semiconductor 
laser and a steering mirror directs the beam. The secondary mirror, the solar 
array and the laser need to be maintained at an acceptable temperature. Hence 
the need for radiators that dissipate the excess of energy that is not converted 
into the laser beam. 

The size of the radiators can be computed considering the steady state ther- 
mal balance between the input power coming from the concentrator and the 
dissipated power through radiation. 

Three radiating areas were considered for the design of the spacecraft: one 
associated to the secondary mirror with area Ar M2 , one associated to the solar 
array with area Ar s , and one associated to the laser with area Ar l . The size of 
each radiating area can be computed from the steady state equilibrium thermal 
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y < 




Figure 11: Illustration of the spacecraft and laser system, showing the two parabolic mirrors 
(Mi, M2), the directional steering mirror (M^), the solar arrays (S) which pump the laser 
(L), and the radiators (R). 



Table 3: Thermal Properties of Spacecraft Elements 



Solar arrays: 




Mirror: 


Radiator: 


Laser: 




Vs a s e s 


T s 


T M2 ct M2 








0.4-0.45 0.8 0.8 


373 K 


373 K 0.01 


0.9 


0.6-0.66 


313 K 



equations: 

Ar s = (a s r]uPiM 2 - - 2e s aA s Tg) j (e R o-T R 4 s ) (32a) 

Ar l = TkVitPiM, (1 - r h ) j (e R( 7T L 4 ) (32b) 

a Rm2 = (a M2 P lM . 2 - 2T^ 2 e s aA M2 ) j (e R crl£ a ) (32c) 

where a s is the absorptivity of the solar array, A$ its area, e s its emissivity, Tg its 
temperature, PiM 2 = t\m Am 1 Sq(t au I ' r 'a) 2 is the input power to the secondary 
mirror, and a is the Stcfan-Boltzmann constant. Then, 775 is the efficiency of 
the solar array, T Rs is the temperature of the radiator associated to the solar 
array, and e R its emissivity. Assuming the efficiency of the laser to be r] L , and 
its temperature T L one can compute the area of the radiator Ar l assuming that 
laser and radiator are in direct contact and that the heat can be transported 
with an efficiency close to 1. This is a reasonable assumption for relatively small 
scale systems that allow the use of a single or bi-phase passive cooling system. 
For large systems a dual phase active system might be required which lowers the 
overall efficiency and increases the system mass. Finally, the secondary mirror 
is assumed to operate at temperature T M2 and has absorptivity a M2 . 

The total mass of the spacecraft is m sc = m^ ry + m„(l + MF t ), where the 
mass of the propellant m p = m^ ry MF v is a fraction of the dry mass m^ry, 
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Table 4: Mass of Spacecraft Elements 



Specific mass: 



g u 0.1 kg/m 2 

g Ud 0.1 kg/m 2 

g L 0.005 kg/W 
£ S 1 kg/m 2 
q k 1.4 kg/m 2 



Mass: 



m& us 500 kg 



Mass fractions: MF H 
MF p 
MF t 



0.2 
0.3 
0.1 



augmented by the mass fraction MF t = 10% to include the mass of the tanks. 

The dry mass rridry = 1.2(m H + m s + m M + m h + m R + mb us ) is the sum of the 
mass of the laser m h , mass of the bus mbus > mass of the mirrors m M , mass of the 
solar array m s , mass of the radiators m R and mass of the harness m H . Given 
the low maturity of the technology employed for this system, we considered a 
system margin of 20% on the dry mass. 

The mass of the harness m H is a fraction of the combined mass of the laser and 
solar array m H = MF H (m s + m L ). The mass of the solar array is m s = lA5g s A s 
where we considered a 15% margin given the high efficiency of the cells. 

The mass of the laser is m L = 1.5g L P L r] L where the margin is now 50% given 
that a semiconductor laser of this size for space applications has not flown yet. 
The mass of the power management and distribution unit dedicated to the laser 
system is included in the mass of the harness while the mass of the bus is 
assumed to account also for the power electronics. The power input to the laser 
is, 



and is a function of the input light power on the solar array, the efficiency of 
the solar array i] s and the reflectivity of the mirrors rj u = 0.90. The loss due to 
power regulation and transmission was considered to be 15% of the generated 
power. 

The mass of the radiators m R = 1.2(A Rs + A RM2 + A RL )g R from Eq. (32) is 
proportional to the area and is augmented by a 20% margin. The total mass of 
the mirror is m M = 1.25(g Md Ad + g M A Ml + g u A U2 ), where A Ml and A M2 are the 
areas of the primary and secondary mirrors. The total mass of the mirrors is 
augmented by a 25% margin given the technology readiness level of the primary 
mirror. 

The thermal properties of the system are reported in Table 3 while the 
values of the specific masses g, mass factors MF and mass of the bus mbus 
arc reported in Table 4. The margins on the mirrors and power system are 
considered to include the marginal use of power to control the spacecraft in 




(33) 
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Table 5: Boundary values on the design variables for the formation design. 



Design Parameter 


Lower Bound 


Upper Bound 


Mirror aperture diameter, d M (m) 


2 


20 


Number of spacecraft, n sc 


1 


10 


Concentration ratio, C r 


1000 


5000 



proximity of the asteroid. As it will be shown later, the required thrust level is 
small and therefore the power demand is marginal compared to the one required 
for the sublimation. 

4-1- Multiobjective Design 

Once the deflection and the spacecraft models are defined, the interest is 
to optimise the formation in order to obtain the maximum value of the impact 
parameter for the minimum mass into orbit, given a warning time. The problem 
can be formulated as follows: 

min (-6) (34a) 
min (n sc m sc ) (34b) 

xG-D 

where the design vector x is defined by [d M , n sc , C r ] T and the search space D 
is defined in Table 5. 

The problem has two objectives, and a mix of integer and real variables. 
MACS was used to solve Eq. (34). The achievable deflection depends on the 
contamination of the optics, therefore the optimisation was run for both the 
shaped and the natural orbits. The result for the case of natural formation 
orbits can be seen in Fig. 12, where the impact parameter is represented against 
the mass of the system and the aperture diameter of the primary mirror for 
a laser with r/L = 0.6 and solar cells with rjs = 0.4, and Fig. 13, where the 
impact parameter is represented against the mass of the system and the aperture 
diameter of the primary mirror for a laser with t]l = 0.66 and solar cells with 
rjs = 0.45. Analogous solutions for the case of the shaped orbits can be found 
in Figs. 14 and 15. 

It is interesting to note that the number of spacecraft increases when the 
aperture diameter increases. This is due to that fact that as the diameter of the 
primary mirror increases the radiator and laser mass increases up to the point 
at which the mass of a single spacecraft exceeds the total mass of two or more 
spacecraft of smaller size. This is a very important point that is in favour of 
the use of a formation instead of a single large spacecraft. 

Furthermore, it has to be noted that the assumption is that the system for 
each spacecraft is scalable. This is actually not true in general as the technology 
for radiators and concentrators cannot be arbitrarily scaled up. In other words, 
technological solutions for small size spacecraft cannot be applied to large size 
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spacecraft without modifications. This is a further reason in favour of the use 
of multiple spacecraft of small size. 

Figures 16 and 17 show the achievable impact parameter for the case of 
the natural formation orbits with two alternative design solutions, a 5 m in 
diameter reflector and a 10m in diameter reflector both with a concentration 
ratio of 5000, i.e., the ratio between the area of the concentrator and the area 
of the spot on the surface of the asteroid is 5000. Figures 18 and 19 show the 
achievable impact parameter for the case of the shaped formation orbits. Figure 
20 shows the sensitivity to the concentration ratio for a fixed warning time of 8 
years. The evident difference between the achievable impact parameter in the 
case of natural and shaped formation depends on the contamination effect that 
stops the sublimation process quite rapidly (less than one year in some cases) 
when the spacecraft rises above the y-z plane. Because the sublimation stops 
at the beginning of the deflection operations, the efficiency of the deflection, in 
the case of the natural formations, is strongly affected by the position along 
the orbit at which the sublimation starts. This is consistent with the results 
presented in Colombo et al. (2009a). 

Laser Efficiency=0.6, Solar Cell eEfficiency=0.4 



1500IK 




Total System Mass [kg] ' Diameter of the Primary [m] 



Figure 12: Total mass of the system against the diameter of the primary mirror of each 
spacecraft and the achieved impact parameter. Natural formation orbits: r\i, = 0.60, r]s = 
0.40. 



5. Effect of Eccentricity 

One may argue that the method is effective only on asteroids relatively close 
to the Sun as the solar collectors need to power the laser. Indeed if the asteroid 
has an aphelion far from the Sun the power can drop below the minimum re- 
quired to sublimate the surface. Using the idea of the shaped orbits, one can try 
to apply the laser concept to asteroids with an increasing aphelion from 1 AU to 
2 AU and with a decreasing perihelion from 1 AU to 0.5 AU. The assumption is 
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Figure 13: Total mass of the system against the diameter of the primary mirror of each 
spacecraft and the achieved impact parameter. Natural formation orbits: r)L = 0.66, r]s = 
0.45. 



Laser efficiency=0.66, Solar cell efficiency=0.45 
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Figure 14: Total mass of the system against the diameter of the primary mirror of each 
spacecraft and the achieved impact parameter. Shaped formation orbits: til = 0.66, r]g = 
0.45. 



that the Earth is moving on a circular planar orbit and the asteroid on a planar 
elliptic orbit. The impact parameter is computed at one of the two intersections 
with the orbit of the Earth and the deflection action starts at the perihelion of 
the orbit of the asteroid. 

Figure 21 shows the achievable impact parameter as a function of radius of 
the aphelion and perihelion for 9 years of warning time, C r = 5000 and a 20 



2G 



Laser Efficiency= 0.6. Solar Cell Efficiency= 0.4 
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Figure 15: Total mass of the system against the diameter of the primary mirror of each 
spacecraft and the achieved impact parameter. Shaped formation orbits: til = 0.60, r]g = 
0.40. 



Mirror diameter d=5m, C r =5000, laser efficieircy=0.6, cell efficiency=0.4 
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Figure 16: Impact parameter as a function of the number of spacecraft and warning time: 
5 m aperture diameter and a concentration ratio of C' r = 5000. Natural formation orbits: 
Vl = 0.60, ris = 0.40. 



m diameter collector. For comparison with the case of Apophis one can notice 
that the achievable impact parameter is substantially high for highly elliptical 
asteroids. There are two good reasons for that. One is that the thrust is applied 
mainly at the pericentre of the orbit but for highly elliptical orbits a variation 
of velocity at the pericentre produces a much higher change of the semi-major 
axis than for low eccentric orbits. The other is that the orbit of the asteroid 
has a much steeper intersection with the Earth's orbit and therefore a small 
variation of the arrival time generates a greater impact parameter. 
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Mirror diameter d=10m, C r =5000, laser efficieney=0.6, cell efficiency=0.4 




Warning Time lyearj 



Figure 17: Impact parameter as a function of the number of spacecraft and warning time: 
10 m aperture diameter and a concentration ratio of CV = 5000. Natural formation orbits: 
Vl = 0.60, r? s = 0.40. 




Figure 18: Impact parameter as a function of the number of spacecraft and warning time: 
5 m aperture diameter and a concentration ratio of CV = 5000. Shaped formation orbits: 
t) L = 0.60, r? s = 0.40. 



If one sticks to the hypothesis used above for the contamination, even in the 
case of natural orbits the spacecraft will experience no contamination as they 
fly above the plume when the sublimation is minimal or null. Therefore, the 
laser ablation seems to be effective even for high elliptical asteroids with high 
aphelion. 
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Figure 19: Impact parameter as a function of the number of spacecraft and warning time: 
10 m aperture diameter and a concentration ratio of CV = 5000. Shaped formation orbits: 
VL = 0.60, vs = 0.40. 




Figure 20: Impact parameter as a function of the number of spacecraft and the concentration 
ratio: 10 m aperture diameter and a warning time of 8 years. Shaped formation orbits: 
r] L = 0.60, T) S = 0.40. 



6. Conclusion 

This paper presented the multidisciplinary design of a formation of spacecraft 
equipped with solar pumped laser for the deflection of asteroids. 

The paper demonstrated that the use of multiple spacecraft is an optimal 
solution to maximise the deflection while minimizing the mass of the overall 
system. In fact as the diameter of the primary mirror increases the radiator 
and laser mass increases up to the point at which the mass of a single spacecraft 



29 



x 10 4 




i p [AU] 



Figure 21: Impact parameter as a function of the radius of perihelion rp and aphelion of 
the orbit of the asteroid. 



exceeds the total mass of two or more spacecraft of smaller size. This is a very 
important point that is in favour of the use of a formation instead of a single 
large spacecraft. A formation, or fractionated system, has the further advantage 
of increasing redundancy and scalability as for a bigger asteroid the solution 
is simply to increase the number of spacecraft. The sizing of the spacecraft 
was based on a simple model in which the mass of the main bus is considered 
constant and the propellant mass is not optimised. These are two limiting 
assumptions that cause an overcstimation of the mass for small systems. At 
the same time the deployment and thermal control systems are assumed to be 
scalable within the range of variability of the design parameters. Looking at 
present technology, this assumption can correspond to an underestimation of 
the mass for large systems. The efficiency of the laser and solar cells are at the 
upper limit of what is currently achievable in a lab environment. Although this is 
an optimistic assumption, current developments are progressing towards those 
limits independently of the deflection of asteroids. It is therefore reasonable 
to expect the system efficiencies presented in this paper in the near future. 
The paper also analyzed the control of the spacecraft in the vicinity of the 
asteroid and showed that with minimal control and propellant consumption the 
spacecraft can be maintained in their desired formation orbits. 

Finally it was demonstrated that the laser ablation concept based on solar 
power is applicable also to high eccentric orbits (deep crossers) with even better 
performance with respect to the shallow crosser case. In fact, for deep crossers 
the deflection action is maximal where most effective, i.e., around the perihelion, 
and the steep intersection between orbit of the Earth and orbit of the asteroid 
amplifies the deflection effect. 
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